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Motivated by recent experiments, we investigate the electron-vibron coupling in suspended carbon 
nanotube quantum dots, starting with the electron-phonon coupling of the underlying graphene 
layer. We show that the coupling strength depends sensitively on the type of vibron and is strongly 
sample dependent. The coupling strength becomes particularly strong when inhomogeneity-induced 
electronic quantum dots are located near regions where the vibronic mode is associated with large 
strain. Specifically, we find that the longitudinal stretching mode and the radial breathing mode 
are coupled via the strong deformation potential, while twist modes couple more weakly via a 
mechanism involving modulation of the electronic hopping amplitudes between carbon sites. A 
special case are bending modes: for symmetry reasons, their coupling is only quadratic in the vibron 
coordinate. Our results can explain recent experiments on suspended carbon nanotube quantum 
dots which exibit vibrational sidebands accompanied by the Franck-Condon blockade with strong 
electron-vibron coupling. 

PACS numbers: 73.63.-b, 73.63.Fg, 63.22.Gh 



I. INTRODUCTION 

The coupling of electronic and mechanical degrees of 
freedom is a promisi ng re search avenue in the physics 
of nanoscale systems .'^'213] 'pj-^g detection of electrome- 
chanical oscillations in suspended graphen^ and carbon 
nanotubes^ represents the current frontier in the field 
of nanoscale resonators. Suspended carbon nanotubes 
(CNT) constitute a particularly interesting model sys- 
tem due to their remarkable electronic properties inher- 
ited from the carbon honeycomb structure, their high 
mobility as well as their many-faceted mechanical prop- 
erties such as low dissipation™ This is highlighted by 
several recent experiments on suspended CNT quantum 
dots which observe pronounced vibrational effects in elec- 
tronic transport.^ 8 9 10 11 12 xhese experiments focus on 
vibrational sidebands in the Coulomb blockade regime 
and suggest the existence of rather strong but sample 
dependent electron-vibron coupling, at least for certain 
vibron modes. Indeed, of the many vibrational modes 
of CNT, only the radial breathing mod^ and the lon- 
gitudinal stretching mode^ESElIlIl have been observed 
in transport experiments. The electron-vibron coupling 
strength seems to vary widely, with some experiments 
showing only weak, if any, vibronic effect^ while oth- 
ers exhibit well-pronounced vibrational sidebands or even 
the Franck-Condon blockade .'^^^ " 

Most of the corresponding theoretical analysis so 
fg^j 8|i3 | i4 tT5] focused on an electron-vibron coupling of an 
extrinsic electrostatic origin. The basic idea behind this 
extrinsic coupling is the mechanical deformation of the 
nanotube induced by the electrostatic interaction be- 
tween the charge on the nanotube and, e.g., a nearby gate 
electrode. These studies leave several issues unresolved 
which we address in this paper: (i) The estimates do not 
succeed in justifying the large values of the coupling con- 
stant A exhibited by some samples, (ii) A criterion to 
determine which vibron modes contribute mostly is still 



missing, (iii) A viable theory must explain the strong 
sample dependence of the electron-vibron coupling. 

In this paper we address these issues by analyzing 
the contribution to the electron-vibron coupling origi- 
nating from the intrinsic electron-phonon coupling in 
CNT. This intrin sic coup ling originates from two prin- 
cipal mechanisms .ISHHiS (j) xhe deformation-induced 
modification of bond lengths induces modulations of the 
electron hopping in the honeycomb lattice, (ii) Distor- 
tions induce a local variation of areas and generate a 
deformation potential. It is interesting to mention that 
these two coupling mechanisms manifest themselves quite 
differently in an effective Dirac equation of the electronic 
degrees of freedom which we employ in this work. While 
the deformation potential enters the Dirac equation as a 
scalar potenti al, the bond length modifications enter as 
a gauge field.' ^^ ' ^^ ' ^^ l 

The coupling constant of the deformation potential is 
estimated to be about an order of magnitude larger than 
the corresponding one for hopping modulations.^ When 
combined with symmetry arguments, this allows one to 
identify which modes couple most strongly to electrons. 
Longitudinal stretching modes as well as radial breathing 
modes involve local area variations and couple strongly 
via the deformation potential, while transverse twist 
modes involve pure shear and thus couple only via the 
weaker hopping-modulation mechanism. Finally, bend- 
ing modes, though coupling via the dominant deforma- 
tion potential, interact only quadratically with electrons, 
since for symmetry reasons their effect cannot depend on 
the sign of the associated deformation. Thus, longitu- 
dinal stretching modes and radial breathing modes are 
the natural candidates to couple strongly to electron s, in 
agreement with experimental conclusions .'^'^Enmnil 

To quantify the strength of the electron-vibron cou- 
pling in suspended CNT quantum dots, we employ 
a mapping to an Anderson-Holstein modeP^I^iElElIll 
which appropriately describes transport through vibrat- 
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ing quantum dots with a discrete electronic spectrum. 
This model assumes that transport occurs through an 
isolated electronic state coupled to few vibrational modes 
of the nanostructure. The corresponding dimensionless 
coupling constant A is given by the shift of the vibron co- 
ordinate (induced by the tunneling electron) measured in 
units of the amplitude of its quantum fluctuations (i.e. 
the oscillator length). An alternative parameter, often 
employed in experimental works is g = X^. To lowest or- 
der, vibron modes will induce vibrational sidebands when 
the corresponding electron-vibron coupling g is of order 
unity. The sideb ands w ill be accompanied by the Franck- 
Condon blockad^^SHH — a low bias suppression of the 
current — when g > I. 

The coupling parameter g does not only depend on 
the speciflc vibron mode but is also sample speciflc, for 
at least two reasons: (i) It is sensitive to the geometrical 
details (radius or length) of the CNT, e.g., through the 
oscillator length of the vibron. (ii) It is sensitive to the 
location of the CNT quantum dot due to the strain proflle 
of the vibron mode. Despite this sample dependence, 
we can give estimates of typical values of g due to the 
intrinsic electron-vibron coupling for the different vibron 
modes of the CNT quantum dots. For the longitudinal 
stretching mode, we flnd upper limit estimates which are 
signiflcantly larger than unity. This is consisten t with 
the observation of vibrational sideband^^EsEHHl and of 
the Franck-Condon blockadd^^l. 

The paper is organized as follows. In Sec. [iT] we 
briefly review several ingredients required to analyze the 
electron-vibron coupling in carbon nanotube quantum 
dots. Specifically, this includes the elastic theory of long- 
wavelength phonons in CNT, the electronic properties of 
CNT in terms of the effective Dirac theory of graphene, 
the formation of CNT quantum dots, and the microscopic 
electron-phonon coupling in CNT. Sec. |III| contains a 
detailed discussion of the dimensionless electron-vibron 
coupling constants of the various vibron modes for CNT 
quantum dots, and contains the central results of the pa- 
per. This section also discusses our results in relation 
to recent experiments. We conclude and summarize in 
Sec. |IV) Some calculational details are relegated to an 
appendix. 



II. PHONONS AND ELECTRONS IN CNT 



A. Vibron modes in suspended CNT 

The acoustic phonon branches of CNT can be de- 
scribed within the elastic theory of a cylindrical mem- 
brane. A CNT with radius r can be viewed as a 
folded graphen e sheet described by the elastic Lagrangian 



with 
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Kstretch — A'^ij + 
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in terms of the mass density poi the Lame coefficients /i 
and A characterizing the in-plane rigidity of the lattice, 
and the bending rigidity k. In Eq. (|2j), u.^ denotes the 
strain tensor of a cylindrical membrane whose compo- 
nents ar^2l! 

h 
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2 u^y = d^Uy 



dyU^ 



(3) 

up to linear order in the tangential and radial distortions 
Mj(r) and h{r), respectively. Here we introduced x and y 
coordinates for the cylinder, fixing the y axis along the 
CNT axis, while the x axis is the curvilinear coordinate 
around its waist (see Fig. [T]) . The strain tensor of a cylin- 
drical membrane contains the additional term h/r in u^^ 
with respect to a flat 2d membrane, describing the strain 
induced by a uniform change in h. In addition, the fl- 
nite curvature of the CNT induces the extra term h/r"^ 
in the bending energy. One readily flnds that the strain 
vanishes for distortions described by 
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u„ ~ sm — 

n \ r 

Uy = constant 



(4) 



h = u cos — J 

with n integer. In particular, the n — 1 case describes 
a rigid translation of the tube in a direction perpendic- 
ular to its axis. The term V^^^^ in Eq. ([T|) describes the 
energy cost to induce a local variation of curvature with 
respect to the cylindrical geometry. This bending energy 
vanishes for the rigid shift Eq. Q with n = 1, but pro- 
duces a gap for the n > 2 modes which guarantees the 
stability of the cylindrical shape of the tube. 

The Euler-Lagrange equations for the Lagrangian Eq. 
([T| and the resulting phononic branches are discussed in 
the Appendix. The distortions involved in stretching, 
breathing, twist and bending modes are depicted in Fig. 
([ij. The corresponding long- wavelength (i.e. qr <C 1) 
dispersions are given by 
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tube. The parameter 7 is thus sample as well as temper- 
ature dependent. Tension modifies the dispersion of the 
bending modes into 
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The dispersion of the n = 1 mode thus becomes linear 
at small momenta and crosses over to a quadratic depen- 
dence at higher momenta. The gapped bending modes 
with n> 2 are less sensitive to tension. 
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FIG. 1: Schematic description of the acoustic phonon modes 
of CNT. a) Unperturbed nanotube. b) Stretching mode, c) 
Breathing mode, d) Twist mode, e) Bending mode with 
n = 1 at finite q. f) Bending mode with n — 2 and q = 0. 



Here, the group velociti es of the stretch ing and twist 
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a^/^/Pq, with A = 2fi + X. The n = 1 bending mode 
has a quadratic dispersion reminiscent of the flexural 
phonon modes in graphene.'^ Unlike the graphene case, 
however, its dispersion is mostly determined by stretch- 
ing energies due to the finite curvature of the nanotube. 
The dispersion of the breathing mode is characterized by 
a finite fre quenc y ujb at long wavelengths. For typical 
parameterSjli^EZl one finds Ar^ ~ 50 eV • [nm] ^ n 
for any realistic nanotube. The bending modulus can 
thus be neglected in the dispersion of the breathing 
mode. The corresponding q = energy gap is then 
Tiuj^ ~ 8 • 10~2 eV/ij_[nm], with Lj_ = 27rr denoting the 
circumference of the tube. As an example, the (10,10) 
armchair CNT with r ~ 7A has huj^ ~ 2 • 10"^ eV, which 
makes the breathing mode frequently too high in energy 
to be excited in low-bias transport measurements. 

In experiments, CNT are clamped by contacts to leads 
at two points separated by a distance L. The wavenum- 
ber quantization q = m'n/ L (with m' an integer number) 
yields discrete phonon oscillator modes, called vibrons. 
These contacts will also frequently induce tension in the 
CNT which can be included into the elastic Lagrangian 
Eq. Q by the modification 
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where 7 measures the tension applied at the ends of the 



B. Electronic states in CNT quantum dots 

Dirac Hamiltonian for CNT. — The electronic proper- 
ties of CNT can be conveniently described in terms of 
the low energy effective Dirac Hamiltonian of electrons 
in graphene^ 



(8) 



with velocity v ~ 10^ m • s ^, where the 2d wavenumber 
k is measured from the relevant Dirac point ±k^, with 
k^, = 27r/(3V3a) (^3, 1) and a = 1.42 1 the bond length. 
The Hamiltonian in Eq. ^ acts on four-component 
spinors (uJi k' k' k' '"a k) Bloch amplitudes in 
the space spanned by the honeycomb sublattices {A/B) 
and Dirac valleys (+/—). The matrices E^, = H^ (g) a^ y 
denote components of a vector S, with Hj and the 
Pauli matrices in the spaces of the Dirac valleys and 
the sublattices, respectively. In CNT the folding of the 
graphene sheet amounts to identifying different lattice 
sites of the same type {A or B) separated by a lattice 
vector of the form T = rfiT-^+fiT^, with = a^S (0, 1) 
and T2 — a-\/3/2 (-\/3, 1) denoting two independent gen- 
erators of the honeycomb lattice.^iS] The corresponding 
nanotube is denoted as [m, n] and it is customary to intro- 
duce the chiral angle = arctan(Ty/T^). Thus armchair 
CNT correspond to 9 = while zig-zag ones stem from 
9 — Tx 12. Owing to the underlying Bloch wavefunction, 
a translation by T modifies the Dirac spinors in the two 
valleys by a multiplicative factor exp[±ik^ • T] . This fact 
can be accomodated in the Dirac description of CNT by 
formally introducing an Aharonov-Bohm vector potential 
a_|_, which induces a fictitious magnetic flux </) given by 
$ = 0/00 = ±k^ • T/27r.li2l Here = / dl • a± with the 
line integral taken along a closed loop around the cylinder 
waist and 0o — h/ei& the flux quantum. Due to flux pe- 
riodicity, we can restrict attention to $ £ [—1/2,1/2]. 
Choosing the a;-axis for the CNT along its waist, i.e. 
X = T, the electronic eigenstates are proportional to 
exp[i(27rZa;/_Lj^ -I- ikyif)], where / is an integer number. 
The low energy Hamiltonian for the CNT thus reads 



hv 



ayky 



(9) 



By direct calculation for a [m, n] nanotube, one obtains 
<i> = ±(m-|-2n)/3. When <!> is integer, the corresponding 
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CNT are metallic characterized by a gapless and linear 
dispersion. Other values of ffi and n lead to semiconduct- 
ing CNT with gapped electronic spectra, the smallest gap 
being given by A = hv^2TT / Lj_. From now on we will 
assume that only the lowest energy subband Z = is 
relevant . 

CNT quantum dots. — The possibility of opening a gap 
in CNT allows for the formation of CNT quantum dots. 
This can be achieved by electrostatic potentials for semi- 
conducting CNT or by inhomogeneities in the graphene 
folding for metallic CNT. In general, the location of the 
quantum dots along the CNT and their length have 
no relation with the overall nanotube size L and can thus 
be sample specific, as indicated in Fig. |2] 
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FIG. 2: Schematic description of the electronic quantum dot 
of size (in gray) along a carbon nanotube of length L. 

In semiconducting CNT electron confinement can be 
induced by an electrostatic potential V{y) (induced in- 
tentionally by gates or by disorder) along the tube direc- 
tion which enters the Dirac Hamiltonian as 



H = 



A + ihvk,, 



A — ifivky 
V{y) 



(10) 



Here we consider potentials smooth on the scale of the 
honeycomb lattice spacing which allows us to restrict at- 
tention to a single cone of the Dirac description. The 
physics of electron confinement can be captured by the 
simplest choice for V{y), namely a square well poten- 
tial of magnitude V{y) = —Vq between y — y^ and 
y — yQ + Lj^ and zero elsewhere, with ?/g the location 
of the left edge of the electronic dot and its exten- 
sion along the tube direction. (Of course, within the sin- 
gle Dirac cone approximation the abrupt variation of the 
potential still needs to be slow on the scale of the bond 
length a.) A potential with Vq > A/2 confines electrons 
in the dot and depletes occupation outside the region 
y G [j/oi2/o + ^d\' yielding quantization of the longitudi- 
nal wavenumber ky = n'vr/L^, with ti' > 1 an integer. In 
the limit hvky <C A one gets a vanishing wavefunction 
for y < yQ and y > y^ + and a non- vanishing spinor 
\n') for y G [yo^Z/o + ij 



1 
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with energy 
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In metallic CNT, quantum dots can be induced by 
a non-uniformity in the folding vector T. If <i> is non- 
integer for y < yQ and y > y^ + L^, and integer for 



y € [j/o'I/o + -^dl' *li^ g^P opening outside the dot con- 
fines electrons within the metallic region with integer <i>. 
The wavefunction for y G [2/0 ' % + then 



^n'(2/) = {VW) 



1 
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in the limit fivky <C A, with a corresponding energy 

n'TT 



E„, = hv- 



L 
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The potential distribution and energy profiles involved in 
the creation of quantum dots in CNT are sketched in Fig. 
|3] While vibron modes of the nanotube involve its entire 



IE 



a) 



kE 



b) 



FIG. 3: Schematic description of the electron confinement 
producing quantum dots in semiconducting (a) and metallic 
(b) carbon nanotubes. Here i/q = 0. 

length, the location and extension of the electronic quan- 
tum dot are essentially determined by inhomogeneities 
and are in principle uncorrelated with the nanotube size. 
This effect is particularly pronounced in dirty CNT. 



C. Electron-phonon coupling in CNT 

The coupling between electrons and phonons in 
graphene and in nanotu bes can be induced by two differ- 
ent mechanisms 

(i) The deformation potential^ is induced by the local 
area variation of the membrane which is associated with 
a phonon excitation. Within the Dirac representation in 
a single v alley, the deformation potential can be written 

j^ jl6ll7l26l 
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(dof) ^ 
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yyj 



(15) 



The corresponding coupling strength has been estimated 
to be as large as ~ 30 eV.^'' 

(ii) Phonon modes involving pure shear do not induce 
a local variation of area and the associated strain has a 
vanishing trace. These modes couple to electrons only 
via the distortion-induced modification of bond lengths 
in the honeycomb lattice. The corresponding modulation 
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in the electron hopping translates in the Dirac language 
into a purely ofF-diagonal term, correcting the electronic 
momentum in the s ame way a s a (fictitious) gauge field, 
with vector potentiaPE^lmlll 



2u 



xy 

- Uyy) 



(16) 



Here t denotes the hopping amplitude, while the rotation 
matrix Z)[— 30] = cos 30 1 + i sin 36* ay is due to the rota- 
tion of the folding vector T by the angle 9 with respect to 
the x — y axes in graphene. The form of the vector poten- 
tial is essentially fixed by the honeycomb lattice symme- 
try. Combined with the diagonal deformation potential, 
the fictitious gauge field yields the total electron-phonon 
coupling matrix 



e—ph 



526 



(17) 



with /[Wy] = 2uxy + i{uxx — Uyy) and the hopping- 
modulation coupling constant 32 = |f — l-5eV.I^ 
The latter is estimated to be much smaller than the cor- 
responding one for the deformation potential. 

Thus, stretching modes and breathing modes induce lo- 
cal variations of areas, yielding a strong coupling via 
the deformation potential, while twist modes involve pure 
shear and couple to electrons only via the weaker " gauge- 
field coupling" . For reasons of symmetry, bending modes 
are qualitatively different, as their effect on the electronic 
energy cannot depend on the sign of the associated defor- 
mation. Indeed, we will show that they yield a vanishing 
expectation value of the strain on the electronic states 
in first order in the displacement fields. Hence, bend- 
ing modes require an analysis of the strain tensor of the 
cylindrical membrane to quadratic order. A lengthy but 
straightforward calculation yields 
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{dxh 
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dxUxf + {dxUyf 



2u., 



dxUy 



dyU^ 



+ {d^h){dyh)-^{dyh)+-{dyU^) + 

+{dxUx){dyU^) + {d^Uy){dyUy) . 

One readily checks that this result respects the require- 
ment that the rigid translation described by Eq. Q with 
n = 1 is strain-free. As will be discussed in detail in the 
next session, for the n — I bending mode a finite cou- 
pling mediated predominantly by the deformation poten- 
tial can be obtained only at finite q, while modes with 
n > 2 couple also for q — 0. Having established the 
form of the microscopic electron-phonon interaction in 
the language of the Dirac hamiltonian, we are now ready 
to systematically analyze the strength of the coupling for 
different types of vibrons in CNT quantum dots. 



III. ELECTRON- VIBRON COUPLING IN CNT 
QUANTUM DOTS 

Electronic transport through quantum dots located in 
suspended carbon nanotubes is adequately described by 
an (extended) Anderson-Holstein model. Motivated by 
experiments, we will focus on the limit in which the vi- 
brational frequencies are small compared to the electronic 
level spacing.l22lThen, we can effectively restrict attention 
to a single electronic state coupled to several vibrational 
modes. We emphasize, however, that the parameters of 
the model can depend on the particular electronic state 
under consideration. 

We now use the results reviewed in the previous section 
in order to evaluate the electron-vibron coupling enter- 
ing into the Anderson-Holstein Hamiltonian for a CNT 
quantum dotPimmill 



leads 



Here 



^vih describes the electronic and vibronic 



properties of the CNT quantum dot, while H^^^^^ + Hr^ 
describes the leads and the tunneling of electrons between 
the leads and the dot. In Eq. (18), is the electron 



number operator with spin s, e the energy of the electron 
state, and U the repulsion energy for double occupancy of 
the dot. The operator b'^q ^ annihilates the vibron mode of 
type u {v = Stretch, Twist, Breath, Bend) and wavenum- 

ber q, with energy hu!q'^\ The electron-vibron interaction 
in can be linear (j = 1) or quadratic {j = 2) in the 
vibronic operators, depending on the vibron mode under 
consideration. 

In order to compute the dimensionless electron-vibron 
coupling constant Xq'^\ we rewrite the electronic Hamil- 
tonian Hc\ in terms of the vibronic displacement operator 



bi-^^)&lq/V2^S 



/2v^ 



f'OSC,q 



\/n/Muj. 



(19) 



and M 



EEt^+^rM- 

s=T,i V 

with the oscillator length 

PqLLj_ the oscillator mass. We can now identify the 
electron-vibron coupling constant Ag"^ with the shift in 
the energy of the localized electronic state induced by a 
static vibron displacement Uq'^^ = lisc.q/V2, measured in 
units of hujq'^^ . This shift can be readily computed in per- 
turbation theory, employing the intrinsic electron-vibron 
coupling Vg_p^ in Eq. (17 1. We find that first-order per- 
turbation theory suffices. Thus, an electronic state \n') 



of the dot is shifted by Ae^ „,(iig ) = ('^'iK-^pftl"' 
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for a given vibron displacement Uq^\ E.g., for a semi- 
conducting dot, the expectation value is taken with the 
electronic state Eq. (111. As a result, we obtain the ex- 
pression 



huj, 



(20) 



for the dimensionless coupling constant 



q.n' ' 

In order to analyze the value of Xg'^^, for the various 
vibron modes, we assume that the nanotube extends be- 
tween y — and y — L, while the electronic wavefunction 
in the dot exists in the interval y G [z/qj 2/o + ^dl' 
picted in Fig. The various vibrons have a quantized 
wavenumber q — to'tt/L along y, with to' an integer num- 
ber (the index q will be replaced by the corresponding to' 
from now on). We start considering the n = modes. 

Stretching mode: The stretching mode implies local 
area variations inducing a deformation potential together 
with hopping modulations which amount, in the Dirac 
language, to the electron-phonon coupling matrix 
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(21) 



Within the qr ^ 1 regime, the stretching mode is de- 
scribed by Uy — usm.{m!'Ky/ L), = and h ~ 
—iqrX/Auy <C Uy. Its linear coupling to electrons is dom- 
inated by the Uyy component in the deformation potential 
(while we neglect u^^ ~ —X/Auyy due to the small ratio 
X/A ~ 0.1), resulting in the electronic energy shift 



1 lT 
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where we introduced 
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Exploiting Eq. ( 20 1 we deduce 



A, 
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Due to the coupling to electrons mediated by the dom- 
inant deformation potential and to their gapless disper- 
sion, stretching vibrons are natural candidates to produce 
a large electron- vibron coupling. The result in Eq. (24) 



shows that the coupling constant can be strongly sample 
specific, as it crucially depends on the geometric details 
of the nanotube and on the location and size of the CNT 
quantum dot (hidden in „,). 



In order to address this point we first consider the 
regime of very clean CNT, in which the dot size coincides 
with the total nanotube length (i.e. y^ = 0, L^ = L). In 

this case we have Im/''n'^'^°~^^ = ^^m' 2n' ^/^' yielding 



(stretch) e 
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i(Stretch) 
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The interference between vibronic and electronic wave- 
functions thus leads to strong selection rules suppressing 
the contribution of most vibron modes. The detectable 
ones are characterized by a coupling constant 



A, 



(Stretch) 
m' ,7i' , clean 
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-^m' 2n' 



(26) 



Only the long wavelength modes in thin CNT can reach 
the strong coupling regime where vibron signatures in 
transport may emerge. 

As the location and size of the quantum dot is not in 
general related to the overall CNT length, we can con- 
sider the opposite limit of a small dot, <C I/, where 
we obtain 



L 



(27) 



This result stresses once more the strong sample- 
dependence of the electron- vibron coupling strength. If 
the dot localizes in regions involving weak vibron-induced 
strain the coupling constant A can be extremely small. 
In contrast, dots localized around the regions of maximal 



strain in the tube can reach values up to Ae, 
gium'TT / L. This leads to a maximal coupling 



(Stretch) 



m' ,n' ,max 



{u) 



A, 



(Stretch) 
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(28) 



where we used the stretching mode dispersion ^ in the 
oscillator length. We point out that vibron branches with 
a linear dispersion have a coupling which depends only 
on the nanotube circumference, not on its length. Re- 
markably, this coupling constant can be quite large for 
thin nanotubes with typical in the nanometer range. 
For them, the measured coupling constant g — X^ for the 
lowest vibron mode can be o f order 10, in agreement with 
recent experiments. '^'^^'^^^ While we focus on the cou- 
pling of vibrons with localized electronic wavefunctions, 
the deformation potential contribution to the coupling 
between the stretching mode and the bosonic excitations 
of a Luttinger model for a very long CNT has been dis- 
cussed previously.!^ This analysis however resulted in a 
very small coupling constant. 

Breathing mode: The breathing mode has a very weak 
momentum dependence in the qr <^ I regime. It is de- 



scribed by 



0, h = u and by the coupling 
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with a dominant deformation potential component. A 
similar analysis to that for the stretch mode above leads 
to Ae^^™'**''^ (u) = Qiu/r and to the coupling constant 



^(Breath) 



7-10- 



I (Breath) 



(30) 



Contrary to the stretching mode, this coupling constant 
depends only on the total tube length but is otherwise 
independent of the vibronic and electronic wavenumbers 
and is not affected by the localization properties of the 
electronic wavefunction. With this respect, the breathing 
mode dimensionless coupling is a much less sample spe- 
cific quantity. Despite the coupling being mediated by 
the deformation potential, the large breathing mode gap 
significantly suppresses the oscillator length as well as 
the final value of the electron- vibron coupling constant. 

Twist mode: Twist modes are associated with deforma- 
tions given by — usm(qy) and Uy = h = 0. Thus they 
induce a pure shear yielding, for a given displacement, a 
hopping-modulation coupling [see Eq. (17)] 



..(Twist) _ / e'^^d 



(31) 



to the electronic Dirac equation. Focusing on semicon- 
ducting dots, this implies 



(Twist), s „ m'-r: ,„„,^m',n' 



As a result, using Eq. (pO| we deduce 



L,, 



(Twist) _ V2ff2^irrm'' C0S(361) /„ 



m'n' 



Twist 



(32) 



(33) 



Similarly to stretching modes, the linear dispersion of 
twist phonons yields a coupling constant dependent only 
on the nanotube circumference. However, the smallness 
of the fictitious gauge field coupling with respect to the 
deformation potential will in the end suppress the twist 
mode electron- vibron coupling. In parallel, the latter 
is sensitive to the CNT chirality via the cos(36') fac- 
tor, yielding a maximal effect on metallic armchair CNT 
(with 9 = 0) and a vanishing one for zig-zag CNT (char- 
acterized by 6* = 7r/2). In the regime of a small dot local- 
ized around the maximum shear, using Eq. (27 1 one can 



estimate the upper value of the dimensionless coupling 



^(Twist) _ 92€c.rn'^ COs(30) 



m ,n .max 



V2fm 



(34) 



Twist 



Replacing the typical parameters of graphene-^- in the 
oscillator length for the twist vibron we get 



, (Twist) 
m' ,n' ,max 



0.17 



cos(36i) 
•y/?n'Lj_[m 



(35) 



Bending modes: The displacements involved in the 
bending vibrons with n > 1 are described in detail in the 



Appendix. It is crucial to notice that at finite momen- 
tum q a non-vanishing component Uy is needed in order 
to correctly obtain the dispersion of the n = 1 mode, 
while the gapped n > 2 branches are well described by 
Mj. = —u/nsiii{nx/r)sm{qy), h = ucos(7ia;/r) sin((7?;) 
and = in the long wavelength regime qr ^ I for 
K ^ Ar^. For n = 1 the and Uyy components 
of the strain tensor up to linear order in the vibron 
displacement are proportional to ((/r)^ cos(x/r) sin((7y), 
while u^y oc (qr)'^ sin(a;/r) cos((7y). Thus they yield a 
vanishing expectation value on the x-independent elec- 
tronic states. In complete analogy, for n > 2 the expec- 
tation value of each component of the strain tensor on the 
electronic states vanishes in linear order in the amplitude 
of the distortions. We thus need to consider the strain 
tensor up to quadratic order as discussed in Eq. |18[ The 
expectation value of u.^y vanishes again, while those of 
u^^ and Uyy do not, giving rise to a non- vanishing defor- 



mation potential. From the j = 2 case in Eq. ( 19 ) we get 
the dominant energy corrections for qr <ti I, 



A (Bond,n>2)/ n 



4^2 



uTim! ^ 




"1 




2L , 


>; 


2 


La 






2 


'1 hr,^ 




) 


2 L 



The soft bending modes with n = 1 behave differently 
from those with n > 2. The former have a gapless dis- 
persion and a big oscillator length, while the latter are 
gapped and have much smaller quantum fluctuations. 
The resulting coupling constants are 



A, 



(Bend,n=l) 



9l 



3(2/1 -I- A) L 



A, 



(Bcnd,n>2) 



9l {n 
K 



I) L 



327r2n4 



L 



(36) 



It has to be pointed out that the functions | ± /a^/ „/ /L^ 
are non-negative for every possible value of m', n' , j/g and 
Ljj, implying a positive coupling between electrons and 
bending vibrons. This can be understood in terms of 
the deformation potential, inducing a positive coupling 
to electrons for increasing local areas, which is always 
the case for bending deformations once compared to the 
equilibrium straight nanotube configuration. In parallel, 
we notice that the electron-vibron coupling for bending 
modes is independent of the graphene mass density but 
is sensitive to both the length and the circumference of 
CNT, as well as to the electronic localization. In the 
case of electronic wavefunctions sharply localized near 
the maximum local strain, the electron-vibron couplings 
reach the maximum values 



A 



(Bcnd,n=l) 
ra'n' ,max 



40 



A 



(Bend,n>2) 
m'.n', max 



L[fim\ 



(37) 
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In terms of transport measurements on suspended CNT, 
due to their quadratic coupling bending modes do not 
produce significant Franck-Condon blockade nor vibra- 
tional sidebands at large bias. We point out that the ap- 
parently large coupling constant of the soft n = I mode 
is in reality suppressed by the fact that all experimen- 
tally realized CNT so far have typical values of Lj_ of 
several nanometers. As an example, the (10,10) arm- 
chair CNT is characterized by Lj_ = 4.4 nm, yielding 

Ai'naT'^'""^'' — L[fj,m\/(2m''^). Only the softest bending 
vibrons in extremely long and thin freely suspended CNT 
could ever reach the regime Amax"^'""^' > 1. 

In addition, in experiments the contact-induced ten- 
sion in the CNT stiffens the low energy dispersion of the 
n = 1 mode and reduces its oscillator length and coupling 
strength. In particular, the coupling constant becomes 

/ ^(Bond) \ 2 

■V (Bend,n— 1) / n— l,g \ (Bcnd,n— 1) /oq\ 

\¥^0 = (Bcnd,Tcns) ' \=0 y'^^' 

with the n = 1 dispersion in Eq. Q . 

We point out that if the CNT shows buckling in its 
equilibrium configuration (e.g. via contact-induced com- 
pressions at the two ends), the symmetry argument for 
bending modes does not hold and a linear electron- vibron 
coupling appears. In principle, it is straightforward to 
treat this situation within the approach presented in this 
paper. However, we refrain from giving explicit results 
here because of the highly non-universal nature of the 
effect. 



on the geometric details of the nanotube such as circum- 
ference or length and, except for the radial breathing 
mode, is sensitive to the location of the CNT quantum 
dot along the vibrating nanotube. 

In very clean CNT, where the quantum dot extends 
over the entire length of the vibrating nanotube, we find 
that the electron-vibron coupling is subject to selection 
rules. This significantly suppresses vibronic signatures 
in the conductance which may be responsible for the ab- 
sence of vibrational features in some experiments on sus- 
pended carbon nanotubes.^ 

In this paper, we have focused on the intrinsic electron- 
vibron coupling in ideal carbon nanotubes. Realistic 
samples may also be subject to defects which affect the 
elastic properties of the carbon nanotubes. Such defects 
are expected to be particularly significant when they 
bring the CNT close to an elastic bistability. Another 
issue which remains an interesting topic for future re- 
search concerns consequences of electronic correlations 
on the electron-vibron coupling. 

Experiments on nanoelectromechanical systems are 
characterized by an increasing ability to design and tune 
the sample properties. We expect that the results of the 
present paper will contribute to extend this ability to 
include the electron-vibron interaction strengths which 
may make previously unexplored parameter regimes ac- 
cessible. 
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V. APPENDIX 

The Euler-Lagrange equations for u^, Uy and h out of 
the Lagrangian density Eq. ([ij at the harmonic level are 



PoUy 



r 

Adluy + fidluy + {A- n)d^dyU^ + ^dyh 



p^h = -A 



\ _ ^ dyUy 



(39) 



with A = 2/i + A. Imposing periodic boundary conditions 
along the x-direction, we look for solutions of Eq. (39) of 



the form w^., w , h oc exp[mx/r -I- iqy — iujt] with q the 
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one-dimensional wavenumber along the nanotube axis. 
For n = we get three equations yielding the twisting, 
breathing and stretching eigenmodes 



PoUJ^Uy = Aq^Uy 



iX-h 



(40) 



Twist mode: The first equation in Eq. ( 40 1 yields the 



twist mode, with — u sm{qy) and Uy = h = 0. Its 
dispersion 



(Twist) 



(41) 



coincides with that of the transverse in-plane phonon 
branch of graphene.^^ 



Breathing mode: The third equation in Eq. (40 1 yields 



the breathing mode, with finite h and 



0. Due 



to the CNT curvatu re, at g = this distortion yields 
a finite gap Wg — ^ {Ar'^ + '«)/ (Po'^^)- <C 1 

regime, we extract Uy from the second equation with uj ~ 
and insert it into the third one, to get the correction 



(Breath) _ 



(l + /3(gr)2 



(42) 



^it^^ ^ = 2(tW { AHi+r,) - 2??) and 77 = n/Ar^ < 1. 
Notice that the group velocity of the breathing mode 
never becomes negative in realistic CNT. 



Stretching mode: The second equation in Eq. ( 40 ) 
yields the stretching mode with finite Uy and — 0. Ob- 
taining h from the third equation in the regime uj <C 
and inserting it into the second we obtain 



Bending modes: Further phononic branches at low 
energy originate from n > 1. Generalizing the deforma- 
tions Eq. ^ at finite q, from the equations (39 1 we get 
the long- wavelength "bending" modes for qr ^ 1 



= sin(na::/r) sm{qy) exp[— icji] 

n 

Uy — —uF^{qr) cos{nx/r) cos{qy) exp[—iujt] (44) 
h = u G„{qr) cos{nx /r) sin{qy) exp[—iujt] 

In the K <^ Ar^ regime the eigenmodes are described by 



F,(.«I)..(l_i(^.^) 



(45) 



A A^ 



and by 



1 



Fn>2ix < 1) 



(46) 



K (n^ - 1) 1 A o 



The resulting long wavelength dispersions are 




"ri>2,gr<gl 



(47) 



X 1 + (qr) 



2 2n^A + n^{A - 2A) + {A + 2A) 
2n2(n4 - 1)^ 



(Stretch) _ 
'■'n=0,q ^ 



lA^-X^ 



(43) 



in lowest order in n/Ar^ <^ 1. 



While the n > 2 modes are gapped by the bending energy 
cost, the n = I mode is gapless, with a soft quadratic 
dispersion predominantly determined by the stretching 
energy Lame coefficients /x and A. 
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